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We reexamine the possibility of d-wave superconductivity in the hole-doped two-dimensional 
Hubbard model. In terms of the gauge field description of the spin fluctuations, we show that d- 
wave superconductivity is unstable in the perturbative region with respect to the on-site Coulomb 
repulsion, U. Whereas in the region where d-wave superconductivity is possible, there is a strong 
constraint on the gap of superconductivity. Analysis of the localized spin moments suggests that 
there is another d-wave (d a .2_ !/ 2-wave) pairing due to the short-range antiferromagnetic correlation. 
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Over the past few decades a considerable number of 
studies have been made on the Hubbard model. In spite 
of its simplicity, we can expect a variety of phenom- 
ena: antiferromagnetism, ferromagnetism, and uncon- 
ventional superconductivity. Among others, the possi- 
bility of d-wave superconductivity has been considerably 
investigated since the discovery of the high T c cuprates. 



The theory of d-wave superconductivity in the Hub- 
bard model is based on the antiferromagnetic spin fluc- 
tuations, [jjj The mechanism is analogous to the p-wave 
pairing mechanism due to the ferromagnetic spin fluc- 
tuations in the 3 He system. |2) In principle, the theory 
relys on the perturbative expansion with respect to U, 
the on-site Coulomb repulsion. 

However, the theory has some disadvantages. First of 
all, it is hard to describe the half-filled case in which 
local spin moments are produced at each site and the 
interaction between the spin moments is proportional to 
the inverse of U. Apparently, this dependence of the 
interaction on U cannot be reproduced by the perturba- 
tive expansion in U. In addition, realistic values of U 
are much larger than the carrier hopping amplitude, t. 
Therefore, we must pay particular attention in applying 
the perturbative expansion in U, especially when we deal 
with the effect of the localized spin moments. 

In this paper, we re-examine the possibility of d-wave 
superconductivity in the hole-doped two-dimensional 
Hubbard model without relying on the perturbative ex- 
pansion with respect to U. Starting from the path- 
integral form of the partition function, we introduce the 
Stratonovich-Hubbard field for the value of the localized 
spin moments. Analysis of it suggests that superconduc- 
tivity based on the antiferromagnetic spin fluctuations is 
unstable in the perturbative region. While in the region 
where d-wave superconductivity is possible, we show that 
there is a strong constraint on the value of the supercon- 
ducting gap. Instead of the d-wave Cooper pairing, there 
is another d-wave (d^^ 2 -wave) pairing due to the anti- 
ferromagnetic short-range correlation. This d 2 .2_ y 2-wave 
pairing state is similar to the spinon pairing in the RVB 
theory. pM However, the crucial difference here is that 



we do not need the U(l) gauge symmetry breaking that 
is essential for the slave-particle gauge theory. jj| 

Formulation- The partition function of the Hubbard 
model is written as Z = J2?c2?cexp(— 5), where 



S = 



dr 



+^E 



(1) 



Here r dependence of all fields is implicit and the summa- 
tion j\ is taken over the nearest neighbor sites. Car- 
rier fields are represented in spinor representation: Cj = T 
( Cjf Cjj. ) and Cj = ( Cjf ) . The on-site Coulomb 



interaction term can be rewritten as, U Y]. 

\2 



(W 4 ) Ej [( n 3l + n ii) ~ ( c i°- c j) + 2 ("jT + n 3i)\ > 
where the components of the vector er = (a x , a y , a z ) 
are the Pauli spin matrices. Introducing Stratonovich- 
Hubbard fields for the charge and spin degrees of free- 
dom, @ we obtain Z = f VcDcDflV^cV^expi-S), 
where the action is given by 



S 







dr 



E Cj (d T - /X) Cj - t ^ (ciCj + CjCi 
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E- 



iU 



7 CjCj 
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up to constant. Here the vector O-,- is a unit vector, 4>j 
represents the value of the localized spin moments, and 
do is the unit matrix in spin space . The scalar (j) C j is as- 
sociated with the charge fluctuations. Note that (j>j = 1 
at half-filling. If we consider hole doping, then <p 3 - takes 
4>j < 1. For these amplitudes, we do not consider the pos- 
sibility of inhomogeneous configurations of them because 
such configurations may compete with superconductivity. 
In order to focus on the possibility of superconductivity, 
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we take the uniform values, that is, <j>j = cf> = const, and 
i(f> C j — —(cjO-QCj) = —(1 — 8), with S the doped hole con- 
centration. As a result, we may write the action in the 
following form 



S= dr 
Jo 



i (id) 
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<t> 2 u 



(3N, 



(3) 



where jl = fj, — U(l — S)/2 and N is the number of the 
lattice sites. Note that the last term in the square brack- 
ets has the form of Hund coupling between the localized 
spin moments and the carriers' spins. 

Gauge field description of the spin fluctuations- Now 
we rotate the spin of the fermion at j-site so that it is 
in the direction of Slj. Such rotation can be done by 
performing the following SU(2) transformation: 



= m, 



(4) 
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and Uj = 


' z >n 
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where, Uj 



Here complex variables z,- CT and defined as Q,- 



ZjCTZj 



Thus, we obtain 



5' 



dr 
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The effective action of the boson fields Zj and Zj is ob- 
tained by integrating out fermion fields fj and fj: 



'off 



-Trln 



O r - M 2~Cz ) °ij 



-UjUiUj + U l d T U l 5 i j\ + 



0N, 



(6) 



perturbative expansion in U 0. Formally, we can re- 
produce the effective interaction of the antifcrromagnetic 
spin fluctuation theory as follows. Expanding the action 
(|6|) in powers of U, the quadratic term yields the boson 
propagator. The effective interaction due to the exchange 
of the bosons corresponds to the interaction derived from 
the antiferromagnetic spin-fluctuation theory JlJ up to 
the paramagnon contributions that leads to p-wave pair- 
ing. H Apparently, we can trust this interaction only in 
the perturbative region with respect to U. In addition, 
the interaction exists only when (f> ^ 0, or coupling to 
the bosons is lost. Nevertheless we will show later that 
<f> = in the perturbative region in U. 

For the analysis of the system, we rely on neither the 
cj)U -expansion nor the t-expansion because both of them 
are reliable only in part of the parameter range of U/t 
and S. Alternatively, we study the system by taking the 
continuum limit. Taking such limit is justified when the 
fluctuations are long-ranged. Since the antiferromagnetic 
spin fluctuations may be long-ranged near half-filling, we 
may take the continuum limit. In the continuum limit, 
the action (Eh is reduced to 



S = dr d 2 ri/j(r,T) 



[d T - fl) (To + lAr — er 2 



~(-ia V + A) 2 
2m 
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where the SU(2) gauge field A^ is defined as An — 
Ha=x y,z^%, a a, = -iUd n U and Sa is derived from 



Eq. (|6|) in principle. 

The system governed by the action (Q) is the fermion 
system with the interaction due to the exchange of the 
SU(2) gauge field A^. Here the gauge field A n is as- 
sociated with the spin fluctuations. The z-component 
At = — i(zj<9 M Zj +zid ti Z[) describes the ferromagnetic 
spin fluctuations. Whereas the x-component A x a = 
—i(z^d a zi — zid^z^) describes the antiferromagnetic spin 
fluctuations. The latter relationship is implied from the 
analysis of the CP 1 representation of the antiferromag- 
netic Heisenberg model. || 

In order to consider the antiferromagnetic spin fluc- 
tuations, we focus on the gauge field AZ- In terms of 

the fields ijj± — ± V\l)/\/2 that diagonalize the gauge 
charge, the action is rewritten as 



where £y = t for the nearest neighbor sites and tij — 
otherwise. If we expand Eq. (^) with respect to <f>U , then 
the second order term yields the RKKY interaction be- 
tween the localized spin moments. Q On the other hand, 
the second order of the expansion of Eq. (||) with respect 
to tij yields the antiferromagnetic Heisenberg Hamilto- 
nian. §,§ 

Before going into the analysis of the action (||), we 
discuss how the formulation is related to the antiferro- 
magnetic spin fluctuation theory that is based on the 



S = J dr J d 2 r IJ2 i> s (r, t) [d T - p, + isA x r 
1 (-iV + sA x f 



i> + (r,T)i/j-(r,T)+^_(r,T)4> + {r,T) }+S A *. (8) 
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<f>U 



The action (j|) has the form of fermions coupled with the 
U(l) gauge field. 
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Analysis of the gap equation- Now we study the possi- 
bility of spin singlet superconductivity based on the ac- 
tion (||). In the following we assume that there exists 
an attractive interaction induced by the exchange of the 
gauge field. In the presence of an attractive interaction, 
the mean field Hamiltonian for the spin singlet pairing 
state may have the following form 



Hmf = 



1 
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The gap A k is evaluated from the gap equation: A k = 



m Ek ^kk- 



tanh 



/3(E k ,+^U/2) 



tanh 



0(E k ,-<t>U/2) 



where £ k = ^£ 2 + |A k | 2 . 
14;k' is derived from Eqs. ( 
gauge field A^. However 
form. 

At zero temperature, the gap equation is reduced to 
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In principle, the interaction 
|) and (0) by eliminating the 
we do not need its explicit 
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Note that in Eq. ([To]) the summation in k'-space is taken 
over under the constraint E^/ > cf>U/2. The presence of 
the constraint on spin singlet pairing states is understood 
as follows. One can see that the second term in the braces 
in Eq. (^) is similar to that of the Zeeman energy term 
produced by applying an in-plane magnetic field to the 
system. In fact, such Zeeman energy term is proportional 
to the applied magnetic field times X)j( c ]t c jl + c i 1 C J T ) - 
(Here the direction of the in-plane magnetic field is cho- 
sen along the x-axis.) Apparently, in the limit of the 
large in-plane magnetic field, spin singlet pairing sates 
are unstable. Similarly, in the large <pU limit, spin sin- 
glet states are not stable. Therefore, if spin singlet su- 
perconductivity is stable, then the superconducting gap 
A should satisfy, 



A > 4>U/2 = A c 



(11) 



In order to find the doped hole concentration depen- 
dence of the constraint, we need to evaluate <fr. For 
the calculation, we assume the staggered form for Clj 
as = (— l) 3 e z , because we are concerned with the an- 
tiferromagnetic spin fluctuations. Note that non-zero 
value of c/> does not imply the presence of the antiferro- 
magnetic long-range order but it just implys the presence 
of the antiferromagnetic correlation because there is the 
phase fluctuations, or the effect of the gauge field A*, as 
well as the quantum fluctuations. 

Now we estimate the value of tj> by solving the saddle 
point equations derived from the action obtained from 



Eq. (El) after integrating out Cj 
with respect to <p yields 



and Ci. The variation 



U 



dj 



vv 3 



1, 
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where a — <pU/ (8t) 



2 sin 2 1 



v = fi/(4£), and K(£) = 
is the complete elliptic integral 



of the first kind. Meanwhile the variation with respect 
to fi yields 



(9) 



Vf 2 — a 2 



djK(^l - 7 2 ) = 1-5. 



(13) 



In deriving these equations, we have used v < and 
\v\ > a, which holds for the hole doped case. 

From Eqs. ( |l2| ) and (|l3| ) we find U/t and S dependence 
of <fi and A c . The boundary between <fi ^ and ef> = is 
shown in Fig. [l] by the solid curve. In the <fr — regime, 
there is no attractive interaction due to the absence of the 
antiferromagnetic spin fluctuations as mentioned above. 
Note that the <j> = regime lies in the smaller value 
of U/t. Apparently, this region contains the parameter 
range of U/t and 8 where perturbation in U is justified. 
Therefore, for the states with <f) — 0, perturbation theory 
is justified. Although we cannot say whether perturba- 
tion theory is reliable in the <f> ^ Q regime, states with 
4> =/= are qualitatively different from those with <fi — 
because the former is unstable in the U/t — »• limit. 
(In the language of renormalization group theory, they 
should belong to different fixed points.) Turning to the 
conditions of d-wave superconductivity, the boundaries 
of A c /t = 0.20 and A c /t = 0.04 in Fig. suggests that 
the occurrence of d-wave superconductivity is restricted 
to extremely small parameter region or we need much 
stronger attractive interaction than the RKKY type in- 
teraction. 

Another d-wave pairing- So far we discuss the possi- 
bility of d-wave superconductivity in the carrier system. 
Now we discuss that there is another ci-wave (d x 2_ y 2- 
wave) pairing associated with the localized spin mo- 
ments. 

In the 4> 7^ regime, we may take the form of the an- 
tiferromagnetic Heisenberg Hamiltonian for the action of 
the localized spin moments. In the Hamiltonian formu- 
lation, it is written as H spm — ' ^i' with 
J = At 2 /U . [0J|] Note that the exchange interaction be- 
tween the localized spins is reduced by factor cf) 2 . Since (j> 
is monotonically decreasing function of the doped carrier 
concentration 6, this reduction suggests the relation to 
the Heisenberg antiferromagnet like behavior of spin sus- 
ceptibility observed in the doped high T c cuprates. |l|] 

Now we discuss another ci-wave pairing. In order to 
describe the localized spin moments ilj, we can intro- 
duce fermion creation and annihilation operators, aj and 

^croj. Due to the constraint |f2^| = 1, the 

,t 



j 1 



fermion operators and aj must satisfy, 
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EsW = i- (w) 

o 

Note that this constraint, the fermion system is half- 
filled, is independent of the doping concentration 5. Un- 
der the constraint (|l4|) , the Hamiltonian H sp i n is reduced 
to, up to a constant term, 

H apia = -lfj2 D lj D *> ( 15 ) 

(id) 

where Z)y = a^aji — OixQjj is defined on each bond. 
Taking the mean fields (Dij) and {D\j), we find that the 
d 2 .2_ J/ 2-wave pairing state and the extended s-wave state 
are degenerate. If we introduce a slight hopping term 
for the fermions, then the d x 2 _ y 2 -wave pairing state is 
stabilized M and the gap is of order of <fi 2 J. Although 
the origin of this d x 2_ y 2 pairing is similar to the spinon 
pairing in the RVB theory, that is, the short-range an- 
tiferromagnetic correlation ||, the crucial difference is 
that we do not rely on the U(V) gauge symmetry break- 
ing that is essential for the slave-particle gauge theory. ||] 
In addition, it should be stressed that this pairing state 
does not imply a superconducting state of the fermions 
because of the constraint (|l4|). 

This pairing state provides another d x 2_ y 2 pairing 
which is independent of d x 2_ y 2-w&ve superconductivity. 
The fact that this pairing state originates from the anti- 
ferromagnetic correlation between the localized spin mo- 
ments suggests that it can be associated with the pseu- 
dogap phenomenon observed in the high T c cuprates. If 
we apply the theory to the Cu site degrees of freedom in 
the CuOi plane in the cuprates, then the d x 2_ y 2 pairing 
can be identified with that observed by angle-resolved 
photoemission spectroscopy, p"^ ] Furthermore, there is 
the experiment that indicates the existence of the pseu- 
dogap of d x 2_ y 2 symmetry also in the insulating phase. 
Jl3| Since our d x 2_ y 2 pairing exists also in the insulating 
phase, the experiment supports the relationship between 
experimentally observed pseudogap with d x 2_ y 2 symme- 
try and the d x 2_ y 2 pairing due to the antiferromagnetic 
correlation. 

Summary- To summarize, we have reexamined the pos- 
sibility of ci-wave superconductivity in the hole-doped 



Hubbard model. We have shown that d-wave supercon- 
ductivity is unstable in the perturbative region in U . 
Whereas in the region where d-wave superconductivity 
is possible, there is a strong constraint on the super- 
conducting gap. Instead, there is another d a .2_ 1( 2-wave 
pairing due to the antiferromagnetic spin correlation. 
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FIG. 1. The boundary between <f> ^ and <\> — in the 
U/t — 5 plane. The boundary is given by the solid curve. The 
dashed curve represents A c /t = 0.20 and the dotted curve 
A c /t = 0.04. 
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